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SKSLL  DESIGNS 


1.  INTRODUCTION 

In  many  industrial  experiments  the  goal  is  to  find  the  settings  of  process 
variables  or  factors  x^,  x^,  xj^  that  will  optimize  some  characteristic  y  of  the 
product.  A  widely  used  strategy  for  such  problems  is  to  approximate  the 
relationship  between  y  and  the  process  factors  by  a  low-order  polynomial.  The 
coefficients  of  the  polynomial  are  estimated  from  data  collected  during  N 
experimental  runs  of  the  process;  the  settings  of  the  x's  for  the  N  experimental  runs 
are  given  by  a  response  surface  design.  A  response  surface  design  for  k  factors  is 
written  as  an  Nxk  design  matrix  D.  To  estimate  the  coefficients  of  the  polynomial, 
the  design  matrix  is  expanded  into  an  Nxp  model  matrix  X  that  has  one  column 
for  each  coefficient  of  the  polynomial  model.  The  estimate  b  of  the  coefficient 
vector  is  then  obtained  from  the  least-squares  formula 

b  =  (X-X)-‘Xy,  (1) 

where  y  is  an  Nx  1  vector  containing  the  response  y  for  the  N  runs.  Several  classes 
of  designs  are  available  for  fitting  the  second-order  polynomial 

k  k  k-l  k 

+  E  E  +  «  (2) 

1=1  1=1  1=1  ;=I+1 

over  a  spherical  region.  The  most  popular  designs  for  fitting  the  second-order 
model  (2)  are  the  central  composite  designs  of  Box  and  Wilson  (1951)  and  the 
designs  of  Box  and  Behnken  (1960).  For  k  =  2,  4,  6,  and  8,  the  uniform  shell 
designs  of  Doehlert  (1970)  require  fewer  experimental  runs  than  the  central 
composite  or  Box-Behnken  designs,  but  the  uniform  shell  designs  have  not  been 
widely  used.  Crosier  (1991a)  proposed  generalizing  the  seven- factor  Box-Behnken 
design  to  a  new  class  of  designs.  These  new  designs  have  many  properties 
(including  the  number  of  runs  required)  in  common  with  Doehlert ’s  (1970)  uniform 
shell  designs;  the  two  series  of  designs  will  be  referred  to  collectively  as  the  shell 
designs. 

The  number  of  levels  of  the  factors  required  by  a  design  is  sometimes  cited  as 
the  reason  for  the  lack  of  use  of  designs  with  many  levels.  Although  designs  with 
few  levels  axe  sometimes  required,  there  are  reasons  to  distrust  the  too- many -levels 
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explanation  of  the  unpopulairity  of  some  designs.  First,  the  design  points  in  a 
completely  randomized  design  are  run  in  a  random  order  so  that  the  factor  levels 
must  be  reset  frequently.  Second,  many  experimenters  who  are  not  familiar  with 
statistical  designs  use  one- factor- at -a- time  designs  in  which  each  factor  is  run  at 
many  levels  (such  individuals  frequently  prefer  designs  with  many  levels).  Third, 
there  is  an  inconsistency  in  the  relationship  of  the  popularity  of  a  design  to  the 
number  of  levels  of  its  factors.  Doehlert  and  Klee  (1972)  show  how  to  rotate  the 
uniform  shell  designs  to  minimize  the  number  of  levels  of  the  factors.  Most  of  the 
rotated  uniform  shell  designs  have  no  more  than  five  levels  of  any  factor;  yet  the 
rotated  designs  are  not  as  popular  as  the  central  composite  designs,  which  have  five 
levels  of  every  factor.  One  reason  for  the  lack  of  popularity  of  the  shell  designs 
proposed  to  date  may  be  their  lack  of  symmetry  rather  than  the  number  of  levels  of 
the  factors  that  they  require.  Section  2  discusses  the  difficulty  of  applying 
asymmetric  designs  to  industrial  problems. 

The  symmetry  of  popular  designs  is  due  to  their  use  of  regular  geometric 
figxires.  The  central  composite  designs  have  factorial  points  that  lie  at  the  vertices 
of  a  k- dimensional  cube  and  star  points  that  lie  at  the  centers  of  the  faces  of 
another,  larger  cube.  The  Box-Behnken  designs  use  only  levels  —1,  0,  and  1  for  the 
factors,  so  the  design  points  must  be  centroids  of  (k  -  m)- dimensional  faces  of  the 
k-dimensional  cube,  where  m  is  the  number  of  ±l’s  in  a  design  point. 

A  regular  simplex  is  the  geometric  figure  formed  by  k  + 1  equally  spaced 
points  in  k-dimensional  space;  an  equilateral  triauigle  is  a  two-dimensional  regular 
simplex.  Both  series  of  shell  designs  are  constructed  from  regular  simplexes.  The 
shell  designs  are  used  for  fitting  a  response  surface  to  k  independent  factors  over  a 
spherical  region,  not  for  fitting  a  response  surface  to  q  proportions  over  a  simplex,  as 
in  mixture  experiments.  [For  symmetric  mixture  designs,  see  Crosier  (1991b).] 
Simplexes  do  not  have  a  unique,  natural  orientation  that  gives  the  best  set  of 
coordinates  to  any  simplex- based  design  in  any  number  of  dimensions.  Several 
orientations  of  simplexes  yield  symmetric  design  matrices  for  the  shell  designs,  but 
the  best  orientation  depends  on  both  the  design  and  the  dimensionality.  Section  3 
presents  a  symmetric  orientation  for  the  uniform  shell  designs  of  Doehlert  (1970). 

Spendley,  Hext,  and  Himsworth  (1962)  noted  that  the  noncentral  design 
points  of  the  seven-factor  Box-Behnken  design  are  the  midpoints  of  the  edges  of 
two  simplexes.  Crosier  (1991a)  proposed  using  the  same  construction  for  other 
numbers  of  factors  and  found  symmetric,  three-level  designs  for  k  =  ll  and  k  =  15. 
Because  of  their  relationship  to  the  uniform  shell  designs  and  their  construction 
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from  simplexes,  Crosier  {1991a)  called  these  designs  simplex-shell  designs,  but  I  now 
prefer  to  use  the  adjective  simplieiai  and  call  them  simplicial  shell  designs. 
Section  4  provides  an  alternative  orientation  of  the  simplicial  shell  designs  that 
yields  symmet-’-  designs  for  all  Jt  Section  5  exzunines  the  efficiency  of  the  new 
designs. 

The  Appendix  gives  coordinates  for  the  vertices,  edge  midpoints,  and  face 
centroids  of  simplexes  smd  derives  the  symmetric  orientations  of  the  uniform  shell 
designs  and  the  simplicial  shell  designs. 


2.  PROBLEMS  WITH  ASYMMETRIC  DESIGNS 

Response  surface  designs  are  usually  applied  by  scaling  the  coded  factor 
ranges  to  the  ranges  of  the  experimental  factors.  A  typical  set  of  experimental 
factors  and  their  ranges  ase  blender  speed  (100-300  revolutions  per  minute,  or  rpm), 
mixing  time  (10-20  minutes),  cooking  temperature  (180-200  degrees),  and  cooking 
time  (40-60  minutes).  There  is  no  difficulty  applying  symmetric  designs  such  as  the 
Box-Behnken  designs  or  central  composite  designs  to  problems  such  as  this  one. 
But  consider  the  difficulties  created  by  a  design  with  unequal  factor  ranges,  such  as 
Doehlert’s  (1970)  four-factor  uniform  shell  design  (Table  1). 

In  coded  units,  the  first  factor  covers  the  interval  [—1,1],  the  second  factor 
covers  the  intervaJ  [— .866,.866|,  the  third  factor  covers  the  interval  [— .816,.816|,  and 
the  fourth  factor  covers  the  interval  [— .791,.79l|.  Two  methods  of  scaling  the  coded 
design  factors  to  the  experimentaJ  fawrtors  will  be  considered  here;  both  methods 
have  shortcomings  for  designs  with  unequal  ranges  for  the  coded  factors. 

Method  A:  Coded  design  diameter  to  experimental  factor  ranges.  Because  the 
noncentraJ  design  points  of  Doehlert’s  uniform  shell  design  aue  at  a  raulius  of  1  in 
coded  units,  Doehlert  and  Klee  (1972)  suggested  linearly  transforming  each  column 
of  the  design  matrix  so  that  ±1  in  the  design  corresponds  'to  the  limits  of  the 
region  of  interest  of  the  natural  variable  in  the  experiment  assigned  to  that 
colunm.'  An  interpretation  of  this  statement  is  that  ±1  in  coded  units  would 
correspond  to  the  low  and  high  values  given  by  the  experimenter  for  each  faictor. 
Thus  —1  would  correspond  to  a  blender  speed  of  100  rpm  and  1  would  correspond  to 
a  blender  speed  of  300  rmp.  Likewise,  —1  would  correspond  to  a  mixing  time  of  10 
minutes  amd  1  would  correspond  to  a  mixing  time  of  20  minutes.  As  —.866  is  the 
lowest  coded  value  '  the  second  factor,  the  lowest  vadue  of  mixing  time  in  the 
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Table  1.  Doehlert’s  Four-Factor  Uniform  Shell  Design 


Point 

** 

^4 

1 

1. 

0. 

0. 

0. 

2 

-1. 

0. 

0. 

0. 

3 

.5 

.866 

0. 

0. 

4 

-.5 

-.866 

0. 

0. 

5 

-.5 

.866 

0. 

0. 

6 

.5 

-.866 

0. 

0. 

7 

.5 

.289 

.816 

0. 

8 

-.5 

-.289 

-.816 

0. 

9 

-.5 

.289 

.816 

0. 

10 

.5 

-.289 

-.816 

0. 

11 

0. 

-.577 

.816 

0. 

12 

0. 

.577 

-.816 

0. 

13 

.5 

.289 

.204 

.791 

14 

-.5 

-.289 

-.204 

-.791 

15 

-.5 

.289 

.204 

.791 

16 

.5 

-.289 

-.204 

-.791 

17 

0. 

-.577 

.204 

.791 

18 

0. 

.577 

-.204 

-.791 

19 

0. 

0. 

-.612 

.791 

0. 

0. 

.612 

-.791 

21 

0. 

0. 

0. 

0. 

experiment  would  be  10.67  minutes.  Similarly,  the  highest  value  of  mixing  time 
would  be  19.33  minutes.  For  cooking  temperature,  ±1  in  coded  units  is  assigned  to 
the  range  180-200  degrees,  but  the  design  has  a  range  of  ±.816  in  coded  imits  for 
the  third  factor,  so  that  cooking  temperature  would  be  varied  over  the  range 
181.84-198.16  degrees  in  the  experiment.  The  last  factor,  cooking  time,  would  be 
varied  over  the  range  42.09-57.91  minutes.  Although  this  method  of  applying  the 
design  is  consistent  with  the  theory  of  optimal  designs,  the  resulting  design,  with 
reduced  ranges  for  most  of  the  factors,  is  generally  imacceptable  to  experimenters. 

Method  B:  Coded  factor  ranges  to  experimental  factor  ranges.  In  method  B, 
the  coded  range  for  each  design  factor  is  transformed  to  the  range  for  the  natural 
variable  (experimental  factor)  assigned  to  that  design  column.  Thus,  —.866  for  the 
second  factor  becomes  10  minutes  of  mixing  time,  .866  for  the  second  factor 
becomes  20  minutes  of  mixing  time,  —.816  for  the  third  factor  becomes  180  degrees 
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for  cooking  temperature,  and  so  on.  This  alternative  method  of  applying  the  design 
solves  the  problem  of  giving  the  experimenter  a  design  with  the  specified  ranges, 
but  it  distorts  the  configuration  of  the  design  points  and  changes  the  mathematical 
and  statistical  properties  of  the  design;  it  is  equivalent  to  rescaling  each  column  of 
the  design  matrix  to  have  a  range  of  dbl  In  coded  units.  Because  the  statistical 
properties  of  the  design  with  rescaled  factors  are  unknown,  statisticiaois  have 
generally  avoided  applying  asymmetric  designs  by  Method  B.  Designs  with  unequai 
ranges  in  coded  units  leave  the  practitioner  with  the  choice  of  a  design  with 
imacceptable  ranges  or  a  design  with  unknown  statistical  properties. 


3.  UNIFORM  SHELL  DESIGNS 

Doehlert  (1970)  obtained  the  imiform  shell  designs  from  a  regular  simplex 
that  has  one  vertex  at  the  origin;  subtracting  each  vertex  of  the  simplex  from  every 
other  vertex  of  the  simplex  yields  the  complete  set  of  design  points.  In  the 
orientation  given  by  Doehlert  (1970),  the  uniform  shell  designs  have  a  different 
range  for  each  factor.  For  k  =  3,  7,  11,  ...,  Doehlert  and  Klee  (1972)  found 
symmetric,  three- level  rotations  of  the  uniform  shell  designs.  Doehlert  and  Klee 
(1972)  give  the  3-  and  7-factor  designs  and  Crosier  (1991a)  gives  the  11-  and  15- 
factor  designs. 

In  the  new  orientation  presented  in  this  section,  the  imiform  shell  designs  can 
be  specified  by  a  few  points;  the  complete  set  of  design  points  is  obtained  by 
generating  dl  permutations  of  the  factor  levels  for  each  basic  point.  Therefore  it  is 
clear  that  (a)  for  any  point  in  the  design,  all  permutations  of  the  factor  levels  of 
that  point  are  also  points  in  the  design,  (b)  for  any  point  in  the  design,  other  points 
with  the  cyclic  permutations  of  its  factor  levels  au’e  in  the  design,  and  (c)  the  design 
points  can  be  grouped  into  sets  with  cyclic  permutations  of  the  factor  levels.  With 
these  ideas  in  mind,  examine  the  new  orientation  of  the  four-factor  uniform  shell 
design  in  Table  2.  The  symmetric  rotation  in  Table  2  can  be  applied  without  the 
difficulties  that  the  rotation  in  Table  1  generates.  The  symmetric  rotation  in 
Table  2  has  an  intuitive  apped  based  on  the  equd  treatment  of  the  factors:  each 
factor  has  the  same  set  of  levels  and  the  same  relationship  exists  between  any  pair 
of  factors.  For  i  even,  the  uniform  shell  designs  can  be  divided  into  two  orthogond 
blocks  (as  in  Table  2).  The  same  number  of  center  p>oints  must  be  used  with  each 
block  to  make  the  blocks  orthogond.  The  orthogond  blocking  for  even  ib  is  a  new 
discovery  and  is  helpful  because  even  k  is  the  case  for  which  the  shell  designs 
require  fewer  design  points  than  the  centrd  composite  or  Box-Behnken  designs. 
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Table  2.  Four-Factor  Uniform  Shell  Design 


Point 

Block 

1 

1 

-1 

0 

0 

1 

2 

0 

1 

-1 

0 

1 

3 

0 

0 

1 

-1 

1 

4 

-1 

0 

0 

1 

1 

5 

-1 

1 

0 

0 

2 

6 

0 

-1 

1 

0 

2 

7 

0 

0 

-1 

1 

2 

8 

1 

0 

0 

-1 

2 

9 

1 

0 

~1 

0 

1 

10 

0 

1 

0 

-1 

1 

11 

-1 

0 

1 

0 

2 

12 

0 

-1 

0 

1 

2 

13 

a 

-b 

~b 

-b 

2 

14 

-b 

a 

-b 

-b 

2 

15 

-b 

-b 

a 

-b 

1 

16 

-b 

-b 

-6 

a 

1 

17 

— o 

b 

b 

b 

1 

18 

6 

—a 

b 

b 

1 

19 

b 

b 

—a 

b 

2 

20 

b 

b 

b 

-a 

2 

21 

0 

0 

0 

0 

1 

22 

0 

0 

0 

0 

2 

a  = 

(3-5^/*)/4 

.191 

b  = 

(1+5V2)/4 

.809 

The  uniform  shell  designs  also  have  the  property  that  if  a  point  is  in  the 
design,  so  is  the  negative  of  that  point.  (Note  that,  in  Table  2,  points  11  and  12  are 
the  negatives  of  points  9  and  10.)  Hence  it  is  only  necessary  to  gpve  hzdf  of  the 
design  points  in  a  table.  To  show  the  pattern  for  odd  Table  3  gives  the  five- 
factor  uniform  shell  design.  Note  that  for  odd  k,  there  is  no  cyclic  group  that 
contains  its  own  negatives  (as  do  points  9-12  in  Table  2). 
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Table  3.  Five- Factor  Uniform  Shell  Design 
Add  negatives  and  center  points. 


Point 

*5 

1 

1 

-1  0 

0 

0 

2 

0 

1  -1 

0 

0 

3 

0 

0  1 

-1 

0 

4 

0 

0  0 

1 

-1 

5 

-1 

0  0 

0 

1 

6 

1 

o 

1 

0 

0 

7 

0 

1  0 

-1 

0 

8 

0 

0  1 

0 

-1 

9 

-1 

0  0 

1 

0 

10 

0 

-1  0 

0 

1 

11 

—a 

6  6 

6 

6 

12 

b 

—a  6 

6 

6 

13 

b 

6  —a 

6 

6 

14 

b 

b  b 

—a 

6 

15 

b 

b  b 

b 

—a 

a  =  (4-6‘/2)/5 
6  =  (1-I-6^/2)/5 

^  .310 
Ci  .690 

To  generalize  the  above  construction  method  to  any  k,  the  )b-factor  uniform 
shell  design  is  given  by  all  permutations  of  the  factor  levels  of  the  four  points 
(— 1,1,0,0,...,0),  {—a,b,b,...,b),  {a,—b,—b,...,—b),  and  (0,0,... ,0),  where 
a  =  (A;  —  1  —  (A;  -I- 1)^/^  )/k  and  6  =  (1 4-  (A:  + 1)  )/k.  As  a  check  on  the  values  of 

a  and  b,  notice  that  a  +  6  =  1.  This  rotation  gives  the  uniform  shell  designs  seven 
levels  (0,±a,±6,±l)  with  the  following  exceptions:  ife  =  3,  o  =  0,  6  =  1,  yielding  a 
three-level  design,  and  A;  =  8,  a  =  6  =  .5,  yielding  a  five-level  design.  The 
two-factor  design  is  a  rotation  of  the  well-known  hexagon  design  and  the  three- 
factor  design  is  the  three- level  design  given  by  DeBaxm  (1959),  Box  and  Behnken 
(1960),  and  Doehlert  and  Klee  (1972).  Note  that  the  five  levels  (—1,  —.5,  0,  .5,  and 
1)  of  the  eight- factor  uniform  shell  design  are  equally  spaced. 

For  seven  factors,  the  three- level  rotation  given  by  Doehlert  and  Klee  (1972) 
will  usually  be  the  preferred  orientation.  For  completeness,  I  point  out  that  the 
seven-factor  uniform  shell  design  has  (at  least)  three  symmetric  orientations:  the 
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three- level  design,  the  seven- level  design  given  above  (using  o  and  b),  and  an 
alternative  seven- level  design  that  has  equally  spaced  levels.  When  the  levels  are 
coded  —3,  —2,  —1,  0,  1,  2,  and  3,  the  design  is  given  by  the  eyelie  permutations  of 
the  four  points  (—3,  0,  1,  —1,  1,  0,  2),  (—3,  1,  0,  0,  1,  2,  — l),  (—3,  1,  1,  2,  0,  —1,  0), 
and  (2,  2,  1,  2,  1,  1,  — l),  their  negatives,  and  center  points. 

4.  SIMPLICIAL  SHELL  DESIGNS 

Simplicial  shell  designs  were  obtained  by  Crosier  (1991a)  by  generalizing  the 
seven-factor  Box-Behnken  design.  The  design  points  are  the  origin,  the  midpoints 
of  the  edges  of  a  regular  simplex  centered  on  the  origin,  and  the  negatives  of  the 
edge  midpoints  (which  are  edge  midpoints  of  another  simplex).  The  simplicial  shell 
designs  share  many  properties  with  the  uniform  shell  designs,  including 
N=k{k  +  l)  +  nQ  design  points  (where  no  is  the  number  of  center  points), 
construction  from  simplexes,  and  symmetric  orientations  given  by  the  permutations 
of  the  factor  levels  of  a  few  points.  Unlike  the  uniform  shell  designs,  the  simplicial 
shell  designs  can  be  orthogonally  blocked  for  any  k:  the  midpoints  of  the  edges  of 
one  simplex  and  center  points  form  one  block  and  the  midpoints  of  the  edges  of  the 
other  simplex  and  center  points  form  the  second  block.  The  same  number  of  center 
points  must  be  used  in  each  block. 

The  midpoints  of  the  edges  of  one  simplex  eire  given  by  all  permutations  of 
the  factor  levels  of  the  two  points  (— 1,-1, 9,9,..., qf)  and  (— r, 6, 6, ...,6),  where 
qf  =  2/(A  — l-i-(A:-t-l)  ^/^),  r  =  (A  — l)/(fc  — H-(A:-i-l)  and  6  is  as  defined  for 
the  imiform  shell  designs:  6  =  (l -f  (ik+ 1)^/^  )//t  Hence  these  points  plus  center 
points,  (0,0,..., 0),  form  one  block;  the  second  block  consists  of  the  negatives  of  the 
points  in  the  first  block.  In  this  orientation,  the  simplicial  shell  designs  for  k  >  3  are 
nine- level  designs.  The  two- factor  design  is  a  seven- level  design  because  the  levels 
±9  do  not  appear;  the  two- factor  design  differs  from  the  uniform  shell  design  only 
by  a  sign  change  of  one  factor.  Table  4  gives  the  two-factor  shell  designs  scaled  to 
cover  the  range  (—1,1).  The  two-factor  uniform  shell  and  simplicial  shell  designs 
are  rotations  of  each  other  and  therefore  have  the  same  orthogonal  blocking 
pattern,  whereas,  for  other  niunbers  of  factors,  the  uniform  shell  and  simplicial  shell 
designs  are  distinct  and  have  different  blocking  patterns.  The  three-factor 
simplicial  shell  design  has  five  levels  but  is  singular. 

For  A:=3,  7,  11,  ...,  symmetric,  three-level  rotations  of  the  simplicial  shell 
designs  exist,  but  (a)  the  three-factor  design  is  singuleir,  (b)  the  seven-factor  design 
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Table  4.  Two-Factor  Shell  Designs 


Uniform  Shell 

SimpliciaJ  Shell 

Block 

Block 

1. 

.268 

1 

1. 

-.268 

1 

-.268 

-1. 

1 

-.268 

1. 

1 

-.732 

.732 

1 

-.732 

-.732 

1 

0. 

0. 

1 

0. 

0. 

1 

-1. 

-.268 

2 

-1. 

.268 

2 

.268 

1. 

2 

.268 

-1. 

2 

.732 

-.732 

2 

732 

.732 

2 

0. 

0. 

2 

0. 

0. 

2 

is  the  Box-Behnken  design,  and  (c)  the  11-  and  15- factor  designs  were  given  by 
Crosier  (ld91a). 

Because  of  the  similar  constructions  of  the  imiform  shell  and  simplicial  shell 
designs,  it  is  possible  to  give  both  shell  designs  for  a  specific  k  in  one  table  by 
redefining  the  constants  in  the  table.  Table  5  gives  the  shell  designs  for  six  factors. 
The  points  listed  in  Table  5  versus  their  negatives  are  an  orthogonal  blocking  of  the 
simplicial  shell  design,  but  for  the  uniform  shell  design  it  is  necessary  to  change  the 
signs  in  rows  16,  17,  and  18  before  assigning  the  points  in  Table  5  to  one  block  and 
their  negatives  to  the  other  block. 

Like  the  uniform  shell  design,  the  seven- factor  simplicial  shell  design  has  three 
(known)  symmetric  rotations:  the  three- level  design  given  by  Box  and  Behnken 
(1960),  the  nine- level  design  given  by  the  construction  method  above,  and  a 
rotation  with  seven  equally  spaced  levels.  One  block  of  the  seven- level  design 
consists  of  center  points  and  the  cyclic  permutations  of  the  four  points  (—3,  0,  0,  — 1, 
0,  -1,  -1),  (2,  1,  1,  1,  0,  -2,  -1),  (2,  0,  1,  -2,  1,  -1,  1),  and  (2,  1,  0,  -1,  1,  1,  -2). 
The  points  in  the  second  block  are  the  negatives  of  the  points  in  the  first  block. 
Designs  with  many  levels  have  been  recommended  for  computer  experiments 
(Welch  et  al  1992),  so  the  nine- level  rotations  of  the  simplicial  shell  designs  might 
be  useful  for  that  purpose.  The  eight-factor  simplicial  shell  design  has  the  simple 
levels  0,  ±.2,  ±.5,  ±.7,  and  ±1. 
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Table  5.  Six-Factor  Shell  Designs 
Add  negatives  djid  center  points. 


Point 


NOTE:  The  uniform  shell  design  has  h=l,  «=0, 
r=.392,  and  6=.608;  the  simplicial  shell  design 
h?s  k=—l,  e=.262,  r=.654,  and  6=.608. 


5.  DESIGN  COMPARISONS 

Lucas  (1976)  compared  response  surface  designs,  including  the  uniform  shell 
designs,  by  their  D-  and  G-efficiencies.  The  G-efficiency,  which  is  based  on  the 
maximtim  prediction  variance  within  the  experimental  region,  is  the  more  sensitive 
criterion:  a  design  can  have  a  high  D-efRciency  and  a  low  G-efficiency,  but  not  vice 
versa.  As  a  practical  rule,  a  design  with  a  G-efficiency  of  50%  or  more  is  good 
enough  for  use.  Table  6  gives  the  number  of  noncentral  design  points  and  the 
G-efficiencies  of  the  Box-BehnVen  designs,  the  central  composite  designs,  and  the 
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Table  6.  N—tiq  and  G-Efficiency  (percent)  for  Designs 


k 

Box 

Behnken 

Central 

Composite 

Uniform 

Shell 

Simplicial 

Shell 

N-rto 

G 

0 

c 

1 

G 

0 

e 

1 

G 

N-no 

G 

2 

8 

96.0 

6 

90.0 

6 

90.0 

3 

12 

71.4 

14 

94.6 

12 

71.4 

12 

0.0 

4 

24 

98.9 

24 

98.9 

20 

70.5 

20 

23.1 

5 

40 

90.9 

26 

87.6 

30 

65.6 

30 

51.8 

6 

48 

67.2 

44 

97.0 

42 

64.4 

42 

76.9 

7 

56 

99.3 

78 

84.7 

56 

62.1 

56 

99.3 

8 

80 

99.8 

72 

61.2 

72 

93.1 

NOTE:  The  central  composite  designs  (CCD’s)  for  5-7  factors  use  a  1/2 
fractional  factorial;  the  8-factor  CCD  uses  a  1/4  frcictional  factorial. 


shell  designs.  The  G-efHciencies  in  Table  6  were  obtained  using  Hq  =  2  center 
points.  The  G-efficiency  of  a  design  is  defined  as  p/V^(x)niax>  where  p  is  the 
number  of  parameters  in  the  model  and  V’(x)max  is  the  maximum  value  of 
F(x)  =  Nx'{X'X)~^x  for  any  point  x  in  the  experimental  region.  The  variance  of 
y  at  X  is  V(x)o^fN,  so  the  G-efficiency  of  a  design  compares  the  maximum  value 
of  V{x)  =  N  Var[y{x)]/(r^  within  the  experimental  region  to  its  theoretical 
minimum,  which  is  p.  The  G- efficiencies  of  some  of  the  uniform  shell  designs  in 
Table  6  are  slightly  lower  than  the  values  rep>orted  by  Lucas  (1976)  becaiise  a  more 
thorough  search  was  done  for  the  point  x  at  which  the  maximum  value  of  K(x) 
occiirs.  Based  on  the  G-efficiency  in  Table  6,  the  four-factor  simplicial  shell  design 
should  not  be  used. 

Both  D-  and  G-efficiency  are  measures  of  information  per  point,  which  is 
supposed  to  allow  a  comparison  of  designs  with  different  numbers  of  design  points. 
From  the  experimenter’s  point  of  view,  the  efficiency  of  a  design  is  primarily 
(perhaps  exclusively)  a  function  of  the  niimber  of  design  points.  Thus  the 
mathematical  efficiencies  are  not  practical  measures  of  how  good  a  design  is,  but 
theoretical  measures  that  indicate  how  well  the  points  20*0  arranged  over  the 
experimental  region.  The  shell  designs  provide  symmetric  designs  with  fewer  runs 
than  the  Box-Behnken  or  central  composite  designs  for  k  =  2,4,  6,  and  8.  For  k  =  3, 
the  Box-Behnken  design  is  the  symmetric  orientation  of  the  uniform  shell  design; 
for  k  =  5,  the  central  composite  design  based  on  a  2^~^  fractional  factorial  is  the 
symmetric  design  with  the  fewest  runs;  and  for  k  =  7,  the  Box-Behnken  design  is 
the  base  design  from  which  the  simplicial  shell  designs  were  derived. 
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6.  SUMMARY 


Response  surface  designs  with  equal  ranges  for  the  coded  factors  have  been 
widely  used,  but  designs  with  unequal  ranges  for  the  coded  factors  have  not  been 
used  frequently.  Asymmetric  designs  give  the  consultant  the  choice  of  obtaining  a 
design  whose  factor  ranges  do  not  match  the  ranges  specified  by  the  client,  or  of 
distorting  the  design’s  geometric  configuration  and  statistical  properties  to  obtain 
the  desired  factor  ranges. 

Symmetric  orientations  were  given  for  the  uniform  shell  designs  and  the 
simplicial  shell  designs.  Both  series  of  designs  have  ib(/b  + 1)  -f-  no  design  points, 
and  are  symmetric  designs  that  require  fewer  design  points  than  the  central 
composite  designs  or  the  Box-Behnken  designs  for  i  =  2,  4,  6,  and  8.  An  orthogonal 
blocking  of  the  uniform  shell  designs  was  given  for  even  L  The  symmetric 
orientations  of  the  uniform  shell  designs  are  generally  seven- level  designs,  the 
exceptions  being  k  =  3  (three- levels),  k  =  S  (five- levels),  and  the  sequence  k  =  7,  11, 
15  ...  (for  which  both  three- level  and  seven- level  rotations  exist). 

The  simplicial  shell  designs  are  a  generalization  of  the  seven-factor  Box- 
Behnken  design.  The  two-factor  simplicial  shell  design  differs  from  the  uniform 
shell  design  only  by  a  sign  change  of  one  factor.  The  three-factor  design  is  singular 
and  the  four- factor  design  is  too  inefiicient  for  use.  The  simplicial  shell  designs 
have  an  orthogonal  blocking  for  all  k,  but  the  simplicial  shell  designs  for  ^  >  4  are 
nine-level  designs  except  for  k  =  7,  which  has  three-,  seven-,  and  nine-level 
rotations,  and  the  sequence  A:=ll,  15,  19,  ...,  which  has  both  three-level  and  nine- 
level  rotations. 


APPENDIX:  DERIVATION  OF  THE  DESIGNS 

Coordinates  for  Simplexes.  The  coordinates  of  the  k  +  \  vertices  of  a  k- 
dimensional  simplex  are  expressed  most  conveniently  as  the  rows  of  a 
(k  +  l)x(A:-)- 1)  identity  matrix  I*+i.  The  distance  between  any  two  rows  of  is 
2^!^  ,  thus  showing  that  does  represent  the  vertices  of  a  regular  simplex.  A 
simplex  centered  on  the  origin  may  be  represented  by  (A;-i- 1)1*+!  —  where 
is  a  (fc -h  1)  X  (fc -H 1)  matrix  of  I’s.  The  coordinates  of  the  vertices  of  the 
simplex  are  therefore  all  permutations  of  the  + 1  values  k,  —1,  —1, . . . ,  —1.  After 
rescaling  to  integral  values,  the  midpoints  of  the  edges  of  the  centered  simplex  may 
be  expressed  as  all  permutations  of  {k  —  l,k~  1,  —2,  —2, . . . ,  —2)  and  the  centers  of 
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the  (;t  —  l)-dimensional  faces  (hereafter,  just  faces)  may  be  written  as  all 
permutations  of  (— fc,  1, 1, . . . ,  l).  Notice  that  the  rescaled  face  centroids  are  the 
negatives  of  the  vertices  (before  scaling). 

Thus  designs  based  on  simplexes  are  more  naturally  expressed  by  k  +  1 
linearly  dependent  variables  than  by  k  independent  variables.  To  obtain 
experimental  designs,  rotate  the  design  matrix  of  A;  + 1  lineau'ly  dependent  variables 
so  that  one  linearly  dependent  variable  becomes  a  constant.  The  vairiable  that  has 
the  constant  value  has  no  effect  in  cailculating  the  distances  between  the  points  of 
the  design;  hence  the  other  k  variables  correctly  describe  the  geometric 
configuration  and  may  be  used  as  the  experimental  design.  It  is  well  known  that 
multiplication  by  m  orthogonal  matrix  O  corresponds  to  a  rotation  if  the 
determinaint  of  0  is  equal  to  1,  and  to  a  rotation  followed  by  a  reflection  if 
det(0)=— 1.  To  reduce  the  Jb  +  1  linearly  dependent  variables  to  the  IVxJk  design 
matrix  D,  it  is  necessary  that  the  orthogonal  matrix  contain  a  constant  colunm.  If 
the  constant  column  is  the  last  coiunm  of  6,  this  process  is  written  in  matrix 
notation  as 

D|dl  =  M0,  (A.l) 

where  M  is  the  Nx  [k  + 1)  matrbc  of  the  design  in  linearly  dependent  variables,  1  is 
an  IVxl  matrix  (column  vector)  of  I’s,  and  d  is  a  constant;  typically,  d  =  0  because 
the  linearly  dependent  variables  sum  to  zero. 

Uniform  Shell  Designs.  For  the  uniform  shell  designs,  the  Nx{k  +  i)  matrix 
M  consists  of  all  permutations  of  (— 1,1,0,..., 0)  md  a  row  of  O’s  for  the  center  point 
(Doehlert  and  Klee  1972) .  For  fc  =  3  (mod  4)  —  meaning  the  remainder  is  3  when  k 
is  divided  by  4  —  Doehlert  and  Klee  (1972)  obtained  symmetric,  three- level 
rotations  of  the  imiform  shell  designs  by  using  a  normadized  Hadamard  matrix  for 
the  orthogonal  matrbc  0  in  equation  (A.l).  [A  Hadamard  matrix  is  a  square  matrbc 
of  I’s  and  — I’s  whose  columns  are  orthogonal  —  that  is,  the  model  matrix  of  a 
saturated  two- level  Plackett  and  Burman  (1946)  design;  the  normalization  consists 
of  scaling  the  Hadamard  matrix  so  that  H'H  =  1.] 

The  symmetric  orientation  of  the  uniform  shell  designs  given  in  Section  3  is 
obtained  by  using  the  orthogonal  matrix 
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(A.2) 


cl 

p  = 

—cl  c 

where  1  is  a  ^xl  matrix  (column  vector)  of  I’a,  /=  (H-(fc+ )/fc,  and 
c  —  {k+  ,  for  0  in  equation  (A.l). 

Simplicial  Shell  Designs.  For  simplicial  shell  designs,  the  matrix  M  of  A;+l 
linearly  dependent  variables  consist  of  all  permutations  of 
(ib  —  1,  Jfc  —  1,  —2,  —2, . . . ,  —2),  their  negatives,  and  the  origin  (Crosier  1991a).  The 
orientation  given  by  the  use  of  a  normalized  Hadamard  matrix  H  for  0  in  equation 
(A.l)  is  the  best  rotation  of  these  designs  for  k  =  7  (the  Box>Behnken  design),  11, 
and  15. 

The  symmetric  orientation  of  the  simplicial  shell  designs  given  in  Section  4  is 
obtained  by  using  the  orthogonal  matrix 

h-g^k  cl 

where  jf  =  (l  — (Ji:  +  1)“^'^^  )/A  and  c  =  (fc  +  l)~^^^,  for  0  in  equation  (A.l).  The 
orthogonal  matrix  F  will  also  yield  nine- level  rotations  of  simplicial  shell  designs. 
The  advantage  of  using  Q  rather  than  P  to  obtain  the  simplicial  shell  designs  is  that 
the  diameter/range  (D/R)  ratio  is  larger  when  Q  is  used  to  obtain  the  designs  than 
when  P  is  used  to  obtain  the  designs.  [See  Crosier  (1991a)  for  a  discussion  of  the 
D/R  ratio.]  The  D/R  ratio  of  the  simplicial  shell  designs  for  3<A:<11  ranges  from 
1.48  to  1.51  for  the  orientation  given  by  Q,  but  it  is  less  than  1.1  when  P  is  used  to 
obtain  the  designs.  Similarly,  the  use  of  Q  to  obtain  the  uniform  shell  designs  yields 
seven- level  designs  with  a  small  D/R  ratio.  In  the  orientation  given  by  P,  the 
imiform  shell  designs  have  a  D/R  ratio  of  2^1^. 

The  rotations  of  the  seven-factor  shell  designs  that  have  seven  equally  spaced 
levels  are  obtained  by  xising  an  orthogonal  matrix  C  for  0  in  (A.l).  To  form  the 
8x8  matrix  C,  start  with  the  7x7  matrix  consisting  of  the  cyclic  permutations  of  (l, 
0,  0,  —2,  1,  1,  0).  Then  append  a  row  of  — I’s  to  form  an  8x7  matrix.  Finally, 
append  a  column  of  I’s  to  form  C.  Normalization  of  C  so  that  C'C  =  I  is 
unnecessary.  The  uniform  shell  design  given  by  C  has  a  D/R  ratio  of  4/3  =  1.33 
and  the  simplicial  shell  design  given  by  C  has  a  D/R  ratio  of  12^^^^  /3  =  1.15. 
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Curtailed  Representation.  Because  the  columns  of  M  are  linearly  dependent, 
one  colunm  of  M  (usually  the  last)  can  be  deleted  without  the  loss  of  information. 
Denote  this  curtzuled  representation  by  M^;  is  not  a  correct  design  for  k  factors 

because  the  distances  among  the  rows  of  are  not  the  same  as  the  distances 

among  the  rows  of  M.  But  the  correct  design  matrix  D  can  be  obt^ed  from  M^. 
For  the  uniform  shell  designs,  D=MjfcA,  where  A  =  I*  —  [(1  +  (it  + )/lt] J  k 
and  Jh  is  a  kxk  matrix  of  I’s.  For  the  simplicial  shell  designs,  D  =  Mj^B,  where 
B  =  Ijk  +  [1/(1  +  (fc  + 1)  J*.  It  is  often  easier  to  form  A  and  B  than  the 

orthogonal  matrices  P  and  Q  on  electronic  worksheets  that  have  matrix  capabilities. 
The  curtailed  representation  of  C  is  the  matrbt  formed  by  the  cyclic  permutations 
of  (2,  1,  1,  -1,  2,  2,  1). 
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